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We report detailed measurements of the low temperature magnetic phase diagram of Er2Ti207. 
Heat capacity and time-of-fiight neutron scattering studies of single crystals, subject to magnetic 
fields applied along the crystallographic [110] direction, reveal unconventional low energy states. 
Er3+ magnetic ions reside on a pyrochlore lattice in Er2Ti207, where local XY anisotropy and 
antiferromagnetic interactions give rise to a unique frustrated system. In zero field, the ground 
state exhibits coexisting short and long range order, accompanied by soft collective spin excitations 
previously believed to be absent. The application of finite magnetic fields tunes the ground state 
continuously through a landscape of non-collinear phases, divided by a zero temperature phase 
transition at /io-ffc ~ 1.5 T. The characteristic energy scale for spin fluctuations is seen to vanish 
at the critical point, as expected for a second order quantum phase transition driven by quantum 
fluctuations. 

PACS numbers: 75.25.+Z, 75.40.Gb, 75.40.-s, 78.70.Nx 
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Despite intense study spanning more than a decade, 
highly frustrated magnetism in the rare-earth-metal ti- 
tanate pyrochlores defies a thorough understanding[lj. 
The pyrochlore lattice is a network of corner-sharing 
tetrahedra, and the archetype for geometrical frus- 
tration in three dimensions [4]. The rare-earth ti- 
tanates in particular have a tendency towards novel 
physics, due to a conspiracy of energy scales where 
dipolar interactions and crystal field anisotropy are 
of comparable importance to the exchange energy 
between magnetic moments. Exotic behaviours at- 
tributed to these pyrochlores include spin icefs'l and 
spin liquid^4j] ground states, magnetoelastic fluctuations 
near zero temperature^, emergent magnetic monopole 
excitations [gI, and order- by-disorder [8[. 

Previous studies of Er2Ti207 identify an antiferro- 
magnetic insulator, with a Curie- Weiss temperature of 
Ocw - -22 K0, g. Crystal field analysis yields XY 
anisotropy such that the magnetic moments of the Er'^+ 
ions, with effective magnitude 3/ib[7], are confined to 
planes tangential to the local tetrahedra on which they 
reside. Heat capacity 9j and neutron diffraction^, 
measurements on polycrystalline samples show a contin- 
uous phase transition to long range order at Tc ^ 1.1 
K in zero field. Using the group theoretical language of 
previous work0, [l3|i the observed magnetic order is a 
six-fold degenerate ground state labeled under the ba- 



sis vector -02 of the irreducible representation T5. This 
ground state is predicted for a simple nearest neigh- 
bour exchange model with planar anisotropy [tII, where 
it has been shown to be selected by thermal fluctua- 
tions in an example of what Villain described as order- 
by-disorder ^l|. However, it is surprising that this sim- 
ple model should accurately describe Er2Ti207, where 
dipolar interactions are expected to be an important per- 
turbation and should select a different ground state |12|. 
Therefore, the possibility that this ground state is se- 
lected energetically by a more complicated microscopic 
Hamiltonian, rather than via order-by-disorder, has not 
been ruled out. 

Single crystals of Er2Ti2 07 were grown using the float- 
ing zone technique and a two-mirror image furnace, in 
3 atm. of air and at a growth rate of 5 mm h~^. Heat 
capacities were measured for a 3 mg crystal using the re- 
laxation calorimetry method with a '^He insert. A similar 
7.3 g single crystal was selected for the neutron mea- 
surements, which were performed using the NG4 Disk 
Chopper Spectrometer [l3| at the National Institute of 
Standards and Technology Center for Neutron Research. 

Heat capacity measurements of single crystal 
Er2Ti207 (Fig. 1) in zero apphed magnetic field 
recreate the previously reported anomaly at Tc 1.1 K, 
which signals the second order phase transition to long 
range magnetic order. Under the application of modest 
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FIG. 1: The phase diagram of Er2Ti207 in an applied [110] 
magnetic field is shown. The contour plot was constructed 
from multiple heat capacity curves, on both warming and 
cooling. A sampling of these is shown in the bottom panel. 
Dotted lines on the top panel trace out the location of the 
maximum at each field magnitude. It is only when H < He 
that this boundary line represents a true phase transition. 
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FIG. 2: Elastic scattering is shown within the [H,II,L] plane 
for panels a) - e). Note that only nuclear allowed Bragg po- 
sitions show intensity throughout the phase diagram. Panel 
f) shows longitudinal cuts through (220) at 50 mK. The exis- 
tence of a broad, diffuse component to the scattering in zero 
applied magnetic field is clear, indicating short range mag- 
netic correlations that coexist with long range order. 



magnetic fields, the anomaly remains sharp, but shifts 
to lower temperature, finally disappearing at a critical 
field value iiqHc slightly larger than 1.5 T. For fields of 
1.75 T and above a broad feature resembling a Schottky 
anomaly is seen, at a temperature that depends linearly 
on the appHed field strength. This is suggestive of a 
crossover to a quantum paramagnetic high-field state, 
reminiscent of that of the prototype transverse field 
quantum critical magnet LiHoF^ (l4. 15. l^. 

Knowing the location of the phase boundaries, we now 
focus on the nature of the different phases. To this end, 
high quafity single crystals of Er2Ti2 07 were studied 
by the time-of-flight neutron scattering technique. This 
allows a simultaneous measurement of both the clastic 
and inelastic neutron scattering cross-section within the 
[H,H,L] plane in reciprocal space. Elastic scattering re- 
veals the magnetic structure of the ordered phases, while 
inelastic scattering informs on the elementary spin exci- 
tations of these states. Pre-existing neutron scattering 
measurements 7, 10] identified the zero-applied-field or- 
dering wavevector as the (2,2,0) Bragg reflection. This 
result is confirmed in Fig. 2, where a dramatic en- 
hancement of scattering intensity at (2,2,0) is observed 
upon cooling below the temperature of the heat capac- 



ity anomaly. However, we also note that a second, broad 
component to the scattering is evident. As is clear from 
Fig. 2 f), when subjected to an applied magnetic field, 
the magnetic refiection becomes sharper and the back- 
ground decreases substantially. This is unambiguous ev- 
idence for the existence of a second length scale in the 
zero-field ground state, where short-range magnetic cor- 
relations coexist with long-range magnetic order. The 
width of this broad feature corresponds to spin corre- 
lations of 15 Angstroms (~ 1-2 unit cells). The origin 
of these short-range correlations is not understood, but 
their presence at T=5G niK is consistent with strong 
quantum fluctuations in the ground state. Muon spin 
rotation results support this, as they observe an anoma- 
lously large relaxation rate below Tc 17 1 . 

Within the [H,H,L] plane, we have measured the mag- 
netic scattering intensity of flve distinct Bragg peaks at 
low temperature (T = 50 mK) as a function of applied 
magnetic fleld (Fig. 3c). In weak fields, the scattering is 
dominated by the large (2,2,0) peak associated with the 
zero magnetic fleld ground state (Fig 3a). An anomalous 
increase is observed in the scattering intensity at the or- 
dering wavevector for ^qH = 0.5 T. This effect has been 
ascribed to the creation of a single domain state in pre- 
vious workQ, however we attribute it to the elimination 
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FIG. 3: (a) Schematic diagram of one of six degenerate 
ground states in zero field, (b) A schematic diagram of 
the high-field ground state. The spins are fully polarized 
within XY planes and the state is non-degenerate. Note that 
spins are subdivided into two sets, residing on perpendicular 
chains and coloured red and blue, (c) The measured field de- 
pendence of the integrated magnetic Bragg intensity of Q=0 
peaks at T=50 mK. The nuclear component to the scattering 
is eliminated by subtracting data taken at 2K. (d) Calculated 
magnetic Bragg scattering is shown for a smooth deformation 
of the magnetic structure taking it from a) to b). A canting 
angle of 0° corresponds to a) while 180° corresponds to b). 



of the short-range correlations that we have observed in 
zero field (see Fig. 2). As we approach He from be- 
low, the intensity at (2,2,0) drops off smoothly as long 
range order is destroyed. Simultaneously, new intensity 
develops at the (2,2,2) position, which dominates the 
scattering in high magnetic field. Strong scattering at 
(2,2,2) is consistent with a highly polarized XY state, as 
illustrated in Fig. 3b. The boundary between the (2,2,0) 
and (2,2,2) dominated phases coincides with the quan- 
tum critical point deduced from heat capacity in Fig. 
1. However, the field-induced reorganization of scatter- 
ing intensity is not sharp at the critical point. Rather, 
a smooth evolution is seen for fields within ^ 1 Tcsla 
of hqHc- We also observe an increase in scattering in- 
tensity at the (0,0,2) position for H He, although this 
magnetic Bragg peak is weak at both low and high fields. 

The field evolution of the magnetic scattering is qual- 
itatively understood as a smooth deformation of the 
ground state from configuration A (illustrated in Fig. 
3a) to configuration B (illustrated in Fig. 3b). To vi- 
sualize this, it is useful to subdivide the pyrochlore lat- 
tice into spins residing on two sets of orthogonal chains, 
coloured red and blue in Fig. 3 a) and b). To go from A 
to B, it is necessary to rotate all blue spins by 90° such 



that they polarize along H, while half of the red spins 
must rotate by 180° and half need not move at all. A 
simple picture emerges in which all spins initially rotate 
through 90°, fully polarizing the blue spins while main- 
taining antiferromagnetic correlations between the red 
spins. This leads to a polarized coplanar state, where 
blue spins point along [1,1,0] and red spins point along 
± [1,-1,0]. Next, while the blue spins remain pinned 
along the [1,1,0] direction, the red spins rotate through 
another 90°, maximizing their projections along the ap- 
plied field direction and smoothly increasing the scatter- 
ing at (2,2,2). This picture gives rise to the calculated 
Bragg intensities displayed in Fig. 3d. 

The above scenario describes the evolution of only one 
of the six degenerate low-field states into the unique po- 
larized high field state. A more realistic picture needs 
to take into account coexisting domains of all six A-type 
states subject to a global symmetry-lowering [1,1,0] mag- 
netic field. In this context, the quantum critical point 
at He is most likely related to a field-induced merg- 
ing of ground states, as recently predicted for Gd-based 
pvrochloresflij]. Nevertheless, we expect the qualita- 
tive description of a manifold of smoothly field-deformed 
states to survive more sophisticated analysis, given the 
striking similarities between the observed and calculated 
Bragg intensities. 

Complementary to our diffraction results, we report 
inelastic neutron scattering measurements of the low- 
energy excitations of Er2Ti207 across the phase dia- 
gram. Excitations have been measured throughout the 
[H,H,L] plane, and a subset of these results are reported 
in Fig. 4. In zero applied magnetic field, the onset of 
magnetic order is accompanied by the appearance of well 
defined spin wave excitations. Two modes are easily seen 
near the (2,2,0) position, one being essentially disper- 
sionless at ~ 0.4 meV while the other goes soft at the 
ordering wavevector, identifying it as a Goldstone mode. 
The dispersionless mode was reported previously, along 
with the absence of the soft mode [3]. It has been ar- 
gued that the unusual ground state of Er2Ti2 07 lacks 
soft spin excitations, or that the lowest energy modes 
fail to couple to the neutronQ. Here, we explicitly show 
that gapless, Goldstone spin excitations exist, since we 
observe spin waves which go to zero energy at (2,2,0), 
(1,1,1) and (1,1,3). This unambiguously resolves a long- 
standing point of confusion: unconventional as it may 
be, the zero-field ground state does indeed support soft 
spin excitations. 

Under the application of finite magnetic fields, the ex- 
citation spectrum displays a rich evolution. At hqH = 
0.5 T, the increase in elastic scattering at (2,2,0) is ac- 
companied by the elimination of the soft mode at the 
same position in reciprocal space. Clearly, the applica- 
tion of such weak fields does not simply create a single 
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FIG. 4: The dispersion of spin waves is shown along the 
[2,2,L] hne in reciprocal space, joining the (220) and (222) 
wavevectors, where Bragg scattering characteristic of the low 
and high field low temperature states is observed, a) shows 
data at T=2 K and H=0, while b) - h) show data at T=50 
rtiK and applied field as indicated. A full spin excitation 
softening is observed for /lo-ff ~ 1.5 T. 



does not correspond to a simple disordered paramagnet. 
Certainly, the difference between the excitation spectra 
of the quantum paramagnetic phase (Fig. 4f), and that 
of the thermal paramagnet (Fig. 4a) is undeniable. In 
larger applied fields, the ground state is dominated by 
Zeeman energy. This is evident at 5 T, where two dis- 
persionless modes are seen. The existence of two such 
modes is indicative of the subdivision of the pyrochlore 
lattice discussed earlier, where the Zeeman energy of po- 
larized red and blue spins are not equivalent due to local 
planar anisotropy. 

To conclude, the low temperature phase diagram of 
Er2Ti207 supports a rich variety of elementary excita- 
tions as the ground state is tuned by external applied 
field. The most dramatic variation with occurs around 
^j-qHc ~ 1.5 T, where a seemingly continuous quan- 
tum critical point is observed. The degree to which this 
landscape of low energy states and quantum critical be- 
haviour can be understood will hinge upon a proper en- 
ergetic description of the system, and we hope this work 
will guide theoretical efforts in that pursuit. 
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domain state, but rather it induces a qualitative change 
in the ground state. In this regime, spin excitations are 
still seen to go soft at (1,1,1) and (1,1,3), but not at 
(2,2,0). Apphcation of stronger fields splits the disper- 
sionless feature into two flat modes, the lower of which 
decreases towards zero energy transfer. This mode soft- 
ening is expected at a second order quantum critical 
point, where the energy gap should vanish as a power 
law in \ H — He \ 1J| . Interestingly, a full mode softening 



is avoided in the canonical case of LiHoF4 due to hyper- 
fine coupling which leaves the quantum critical region 



inaccessible [16[. This does not occur in Er2Ti207. At 
jj,oH = 1.5 T, near He, the system is ungapped to within 
our instrumental resolution. An intriguing picture of co- 
existing phases emerges, where critical fluctuations of 
the low field state are observed around (2,2,0), accom- 
panied by a precursor of the intense dispersive excitation 
characteristic of the high field state. This mode is clear 
at 2 T and 3 T, showing that the high field ground state 
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